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Abstract We consider a complex scalar field minimally coupled to gravity 
and to a U(l) gauge symmetry and we construct of a first order symmetric hy- 
I perbolic evolution system for the Einstein-Maxwell-Klein-Gordon system. Our 

analysis is based on a 1-1-3 tetrad formalism which makes use of the components 
of the Weyl tensor as one of the unknowns. In order to ensure the symmetric 
i hyperbolicity of the evolution equations, implied by the Bianchi identity, we 

introduce a tensor of rank 3 corresponding to the covariant derivative of the 
Faraday tensor, and two tensors of rank 2 for the covariant derivative of the 
O ,' vector potential and the scalar field. 
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CN) ' 1 Introduction 

> ■ 

In this article we discuss the construction of suitable evolution equations 
for a self-gravitating charged scalar field governed by the so-called Einstein- 
I Maxwell- Klein- Gordon system. It is well known that General Relativity admits 

initial value problem formulation whereby one prescribes certain initial data on 
a 3-dimensional hypersurface, and one purports to reconstruct the spacetime 
■ associated to this initial data — a so-called Cauchy problem. The formulation 

of an initial value problem is a natural starting point for a wide variety of an- 
alytical studies of the qualitative properties of the solutions to the equations 
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and a necessary starting point for the construction of the numerical solutions. 
Examples of qualitative aspects of the solutions requiring an suitable initial 
value formulation are the discussion of local and global existence problems and 
the analysis of the stability of certain reference solutions. 

In this article we approach the construction of the evolution equations of 
the Einstein-Maxwell-Klein-Gordon system from the point of view of mathe- 
matical Relativity. Hence, the amenability of our analysis to analytic consider- 
ations takes precedence over numerical considerations. Our discussion is chiefly 
concerned with the construction of a first order system of quasilinear hyper- 
bolic evolution equations (FOSH) out of the Einstein-Maxwell-Klein-Gordon 
equations. In general, FOSH systems can be written in the following form: 

A°dtv - A'djw = Bv, (1) 

where the variables of the system are collected in a n-dimensional vector v, t 
denotes a suitable time coordinate and x ~ {x^.x^^x^) denotes some spatial 
coordinates and j = 1, 2, 3. The matrices A'^ and A-' are matrix valued func- 
tions depending on the coordinates (t, x) and the unknown v — that is, the 
system (1) is, in general, quasilinear. The system is said to be symmetric hyper- 
bolic if the matrices A'' and A^ are symmetric and if A*^ is a negative-definite 
matrix. 

As a consequence of the construction presented in this article, one auto- 
matically obtains a local existence and uniqueness result for the equations 
of the Einstein-Maxwell-Klein-Gordon system. The symmetric hyperbolicity 
of the resulting evolution equations ensures that the Cauchy problem for the 
systems is well-posed. In other words, if one prescribes suitable initial data 
on an initial hypersurface, a unique solution exists in a neighborhood of that 
hypersurface — local in time existence. The solutions depend continuously on 
the values of initial data [1]. The question of the well-posedness of the evo- 
lution equations of the Einstein-Maxwell-Klein-Gordon system is a problem 
that touches upon many aspects of current theoretical and numerical analy- 
sis of many physical phenomena. An suitable point of entry for the extensive 
literature in this topic, with particular emphasis on numerics, is given in [2]. 
Scalar fields enter as important ingredients in many theoretical models of con- 
temporary physics. Here, we consider a self gravitating scalar field, minimally 
coupled with the gravitational field via the Einstein equations and the electro- 
magnetic field by the coupling constant q — the so-called scalar boson charge. 
In the presence of strong gravitational fields, scalar fields are described by the 
general-relativistic field equations. These configurations arise in many areas of 
high-energy Cosmology as inflationary models carried by scalar fields. More- 
over, scalar fields are also candidates for dark matter source. In astrophysics, 
scalar fields arise in stellar models as configurations of hypothetical scalar 
particles known as bosons stars, or as a boson core of very compact stars [3]. 

Our analysis is based on a 1 + 3 tetrad formalism. The hyperbolic reduction 
procedure described in the present article borrows from the discussion of the 
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evolution equations for the Einstein-Euler system by H. Friedrich in [4] — see 
also [5,6,7] and the generalization to Einstein- Euler-Maxwell system in [8]. 
In this reference a Lagrangian gauge was used to construct the required hy- 
perbolic reduction and to obtain the desired evolution equations. The central 
equation in this discussion is the Bianchi identity. It provides evolution equa- 
tions for the components of the Weyl tensor. The addition of electromagnetic 
interactions to Friedrich's system in [8] destroys, in principle, the symmet- 
ric hyperbolic nature of the evolution equation as derivatives of the Faraday 
tensor enter into the principal part of the Bianchi evolution equations. This 
difficult was handled by the introduction of a new field unknown corresponding 
to the derivative of the Faraday tensor for which suitable field and evolution 
equations can be obtained. This strategy can be adapted to the case of a com- 
plex scalar field case coupled with the electromagnetic field. In addition to the 
auxiliary variable associated to the covariant derivative of the Faraday tensor, 
our analysis required the use of two further tensors of rank 2 for the covariant 
derivative of the vector potential and the scalar field. The Lagrangian approach 
to the description of a scalar field follows the ideas developed in [9] for a real 
scalar field. It is important to observe that the extra gauge freedom induced 
by the frame representation of the Einstein-Klein-Gordon system used in the 
present article is associated to the evolution of the spatial frame coefficients 
along the flow of the time-like frame — in particular, in order to implement 
the Lagrangian description it is required that the timelike vector of the or- 
thonormal frame follows the matter flow lines. In order to do this the scalar 
field must satisfy suitable regularity conditions. In fact, the local nature of the 
treatment (fixing suitable initial data existence and uniqueness of a solution 
to the evolution equation can be established locally in time) assures the scalar 
gradient to be timelike only in a neighborhood of the initial data. Afterwards 
this conditions may no longer be satisfied and the the scalar gradient can be 
spacelike or null, and in this case the entire set up breaks down. The remaining 
frame components are chosen to be Fermi propagated along this direction. 

The present article is structured as follows: In Section 2 we write and 
discuss the relativistic equations describing a charged scalar field. The tetrad 
formalism used in this article is briefly reviewed in Section 3. General remarks 
concerning the reduction procedure to obtain suitable evolution equations are 
given in Section 4; the resulting evolution equations are discussed in Section 
4.2 and the subsequent sections. The auxiliary fields and the gauge conditions 
are presented in Section 4.6. A summary of the evolution equations is given in 
Section 4.12. Some concluding remarks are given in Section 5. 



2 The charged scalar field equations 

In the present article we will consider the Einstein field equations 
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with matter source given by a self-gravitating charged (i.e. complex) scalar 
field minimally coupled to gravity and to a U{1) gauge field A^. The total 
energy momentum tensor for this system is given by 

where, using units such that h ~ c = 1, 

T^r' = ^g^.F^.F^P - .g'^F^A^.p, (4) 
denotes the energy momentum tensor of the free electromagnetic field and 

T<^<^> = {D^<Pr {D,-P) + {D^<P) {D,<Pr - g^^.g^" {Dx-P)* {Dp<P) (5) 

is the energy momentum tensor for the charged scalar field — see e.g. [10]. In 
the previous equations Z?^ = V^^ + ig^p, where the constant q is the boson 
charge and stands for the Levi-Civita covariant derivative of the metric 
of signature (+, — , — , — ). Moreover, * denotes the operation of complex 
conjugation while 

F^^ = V^A^ - V^A^ (6) 

is the electromagnetic field tensor (Faraday tensor). 

In the sequel, it will be convenient to write the scalar field in terms of 
two real scalar fields 9 and such that 

${x) = e*^(^V(2;), (t^* =0, 0* = 9. (7) 

The energy momentum tensor of equation (5) then reads 
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-g^. (^g^''Vx<l>Vp^ + a^ + q^A^ + 2qaxAA ] , (9) 



where the vector field (t^ is defined by cr,, — V^0, and wc have used the 
notation A^ = A\A^. Consequently, it is possible to write the tensor T^^°' as 

= Tj^) + yx"^^ + U'^nv + U''TiT\ (10) 



where 



r/,^) = V^0V,0 - ^ffp^g^^VAc/'Vp^, (11) 

T/-) =2a^a, -.g^,a2, (12) 

Tll^ =q^2ApA,-gp,A^), (13) 

T/r^ = 2g(2a(^A,) - g^.a^A^). (14) 
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Consistent with the requirement that V^T^ii/ = 0, wc require the scalar field 
^ to satisfy the Klein- Gordon equation 

D^D>'^ = 0. 

Equivalently, one has that 

V^V^-?* + 2iA''V^<P + i.g^'^V^A^'P - A^^ = 0. (15) 

In what follows, for simplicity of the presentation, we set q — 1. 
Written in terms of the real fields (j) and ct^, equation (15) reads 

^V^V^0-aX-^2-2A^a„ = O, (16) 

V^a^ + \c7^,\7^(b + V^V^0 + ^A^V,,^ + 5^"(V^A,) - 0. (17) 



Finally, the electromagnetic field is described by the Maxwell equations in 
the form 

^l^.F,x] = 0, VF^, = 02 (a^ + A^). (18) 
3 Tetrad formalism 

In the present article, the Einstein fields equations will be expressed in terms 

of a frame formalism introduced in [4]. To this end, let {ea}a=o 3 denote a 

basis of frame vectors on the spacetime (A1,5^iy) satisfying gab = g{e.a,Gb) = 
Vab = diag(l, — 1, — 1, — 1). We denote by (w") the corresponding dual basis 
(cobasis). Here, and in the rest of the article, Latin letters a, &, . . . are used as 
spacetime frame indices taking the values 0, . . . , 3 while Greek letters /i, v, . . . 
denote the tensorial character of each object — i.e. they are spacetime indices. 
The Latin letters i, j, k ... will be used as spatial frame indices taking the 
values 1, 2, 3. 

The frame fields Ca and the cobasis are expressed in terms of a local 
coordinate basis as 

ea = e,^a^, = Lj%dx^. (19) 

Thus, 

= u^lel = <5/, (20) 

so that the metric tensor can be written as 

g^. = Vabu^y., 9^" = V'^'e'^el (21) 
The commutation coefficients D^^ for the tetrad e" are defined by 



[ea, Cb] = DJec ^ {e^d^b ~ <dX) e. 



(22) 
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Finally, the connection coejjicients (or Ricci rotation coefficients) 1]^° for the 
tetrad e" are defined by the relations 

Vaeb = r,\e,, Vauj' = -r,\Lo\ (23) 

Hence, equation (22) can be rewritten as 

[ea,e,] = (r/,-n=Je,. (24) 

In particular, since Ca ('76c) = 0, one has that 

ra[bc) = 0. (25) 

The components of the Riemann tensor with respect to the frame Ca are given 
in terms of the connection coefficients by 

The components of the Riemann tensor admit the decomposition 

Rabcd = Cabcd + {ga[cSd]b — 9b[cSd]a} i (27) 

where Sab denotes the Schouten tensor 

Sab = -Rah — ^Rgab, (28) 

with Rab = R'^acb components of the Ricci tensor and R = g°'^Rab the Ricci 
scalar. Finally Cabcd denotes the components of the Weyl tensor with respect 
to Ca- The components of the curvature tensor satisfy the Bianchi identity 

V[ei? Hrf = 0. (29) 
As it is well known, the contracted version of the above identity leads to 

V°Ga6 = 0, Gab ^ Rab — -^Rgab, (30) 

where Gab denotes the components of the Einstein tensor. Defining the Friedrich 
tensor Fabcd via 

Fabcd = Cabcd — ga[cSd]b, (31) 

one concludes, furthermore, that 

Fahc = ^dF'^abc = 0. (32) 

Taking the Hodge dual of equation (32) with respect to the index pair cd, we 
obtain another equation of the form (32) for a tensor Fabcd defined by 

Fabcd = Cabcd + T^^pb^^ acd, (33) 

where Cabcd denotes the components of the completely antisymmetric Levi 
Civita tensor with respect to the frame Ca and C*^,^^ = ^Cabef^'^^cd — 
e.g. [4,11,12,13]. 
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4 The hyperbolic reduction procedure 

4.1 General considerations regarding Friedricli's frame formulation of the 
Einstein field equations 

Following [4], it is convenient to introduce here the following notation: 

N'' = 5^, = iV^e, = eo. (34) 

A tensor Ta^...a is said to be spatial if any contraction with with vanishes. 
The subspaces orthogonal to A^" inherit the metric hab ^ gab — NaNt, (indices 
are raised and lowered using gab)- Thus, ha^ is the orthogonal projector into 
these subspaces. 

For a given tensor, any contraction with N will be denoted by replacing 
the corresponding frame index by N. The projection of a tensor with respect 
to hj' will be indicated by '. Thus, a tensor Tabc one has that 

TaNb = hJ^NPhb'^T^p,. (35) 

A spatial vector satisfies Ta^-.-ap ~ T'^^...^^. In what follows, let tabc = ^Nabc^ 
where 69123 = 1. In terms of the latter one has the decomposition 

eabcd 2 {N[aeb]cd - f^ab[cNd]) ■ (36) 

Given a spatial tensor we define the spatial covariant derivative 

VaTa,...ap = V bTb,. . -b^K^ ha,""' ■ ■ ■ /la/"' (37) 

In particular, it can be readily verified that Vahbc — 'Dabbed — 0. 
In the sequel, it will be convenient to introduce the tensors 

a"-iV''VbiV^ Xab = ha'V,Nb, x = h''\ab. (38) 

These are related to the connection as can be see from 

VaN' = Naa'' + Xa'', a'' = W'r^^, Xab = -ha''hb''r^a- (39) 

For future use, we notice the decomposition of the Weyl tensor as 

Cabcd = 2(/fc[,^d]a - la[cEd]b) - 2(7V[,^rf]pe^'^, + N^aB^pi^J (40) 

in terms of its electric Eab = C'/vaArfc ^-^^ magnetic part Bab ~ C!*'NaNb relative 
to Na, where lab = hab ~ NaNb- The Schoutcn tensor Sab will be expressed in 
terms of the matter fields. Accordingly, we write 

Sab = ST + Sl:-\ (41) 

where 

ST-<Kr~l9abg''V)^ (42) 
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with a similar expression for -S^™'. Using equations (40)and (41) in equation 
(27), the Riemann curvature tensor Rabcd, can be written in the form 

Rated — R^bcd + ^abcd + ^abcd^ (43) 

where 

^abcd = '^({E'^ + -^^)(^ {9a[c9d]b - gb[cgd]a) " 2iV[6.gQ] [c^^d] ) 

-2E[bga][cEd] + 2B[bga][cBd] + 2{N[agb][cGd] - 5[a36][c^d])^ ,(44) 

K!bcd = '^{h[cEd\a - h[cEd]b) ~ '^{^[cBd]pf^ab + ^[aBb\p^^ cd) ^ (45) 
B^bcd = ga[cSd]b ~ 9blcSd]^, (46) 

with defined in equation (42). Finally, from the Bianchi identity for 

the tensor Fabcd, equation (32), we obtain the decomposition 

Fabc = Na [F'^bNNc - F'^.^Nb) - 2Fi^^,N,^ + N^F'^^c + ^a^o (47) 
where 

E'aNb ~ ^NF'jq^j^b + C^i^caJVb ~ {E'Nacb + E^aNb) + ^a-^WAfJVb 

~X'' F'cadb ^ XaENcNb ^ Xb'^^ivaTVc + X'^aE'cNNb + XE'^aNb (48) 

with Cn denoting the Lie derivative in the direction of N"". From equation 
(31) we compute the following components for the Fricdrich tensor Fabcd and 

Fabcd- 

EjvNNa = 0, -P'ivWAfa = ~ -^TecN'^ha'' , (49) 

EwaNb = Bab, E^aNb = Eab - -^Tdchb'^ ha" , (50) 

FNNab = ^{TpvN'")el^bN''', F'j^^ab^O, (51) 

FaNNb = -Bab + | (Tp^iV^) eP aubN^ , (52) 

FaNNb = -Eab + ^hab [T^^iVA^"] , (53) 

FNabc = Eape^bc + ^ {FpvK") E^ubc^", -F'ivabc = -Bap^^bc, (54) 

^aiV6c = -^ape^c + ^ (T.pTV") eP dv f ha" hc^ , (55) 

^aATftc = ^apE^c " | [/^afc {TfuN^hJ) - (T/„iV"/ifc-^)] , (56) 

i?^,^, = -2E^[b^P ^ + ^Tp^ahb"ePaucN\ (57) 

F'abNc — —Bcp^^ab + -^hacTuvN"" h]!, , (58) 

^I^^c. = -^p,e^afce'crf + f (Tp^/^ffc) eL/^"/^c^/i/, (59) 

Fabcd = 2(/h[ci?(j]a - h[cEd]b) - I^Tuvhb'hc^^ahd]' ■ (60) 
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4.2 Evolution equations for the eleetric and magnetie fields 
The Maxwell equations are given by 

^[aFbc] - 0, V'Fab = Jt, (61) 

where 

Ja = -^^{<7a+Aa). (62) 

As it is well known, the electromagnetic field tensor (Faraday tensor) Fab 
can be split into its electric and magnetic parts, Ea = FabN^ and B° = 
i^abcdjy^p^^^ Components with respects to a flow. More precisely, one has that 

Fab = 2E[aNb] - eabcdB'N^, (63) 

where NaN'^ = 1. Using the decomposition into electric and magnetic parts, 
the electromagnetic energy-momentum tensor of equation (4) can be written 
as 

Tir' ^ -\NaNb{E'' + B'^) + '^{E' + B'^) + Pab - 2g^aNb), (64) 

where we have written E'^ = EaE"' and B^ = B'^Ba, and Pab denotes the 
symmetric, trace-free tensor given by 

Pab^^iE^-+B^)-iEaEb + BaBb), (65) 

and 

Ga = eauvdE^B^'N^, (66) 

denotes the Poynting vector. Projecting equations (61) along direction longi- 
tudinal and transverse to the vector A^^,and orthogonal to the N'', one obtains 
the Maxwell evolution equations 

E^f^ = -2E''hf[aVb]N'' - eabcdhf^'V^'iB^N'') ~ h.Kh, (67) 
= -2S"7i;[,Vb]Ar^ + eabcdhf''V''{E'N''), (68) 

and the Maxwell constraint equations 

D'^Ba = -eabcdN^E'^V'^N'', (69) 
D^Ea = eabcdN''B''V^N'' + N^Jb, (70) 

where E = N°-VaEb stands for the covariant time derivative of Ea along 
N, W{a) = ha^'wb is the orthogonal projection of covariant time derivatives 
along the A^*^ and DaWb = h^h^'^uWy denotes the fully orthogonally projected 
covariant derivative of a vector Wa- Finally, we note that 

eabcdhf''V''{X''N'') = -cnAXf + eafcdN'^X'^N'', 

where curlA/ = efacdN'^VB". 
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4.3 Evolution equation for Bat 

The evolution equation for the magnetic part of the Weyl tensor, Bab is en- 
coded in the component Fl^a\N\b) Friedrich tensor. More precisely, one 
has that 

= P{a\N\b) = ^NBab - D dE c(a^b)'^'^ + '^O.ce.'^'^ (aEb)d " (aBb)c 

-2x(a'4)c + xBab - XcdBpqeP^ iae^\) + (71) 
where -F( = -f^(V|f,) + ^HN\b) '^^'^ written as 

+ ^eP,uibX'a)N^Tp,N\ (72) 

Substituting equations (3) and (4) in equation (72) wc find that the electro- 
magnetic contribution to the evolution equation of Bab is given by 

KaWb) ^ -^{eV^{N''eP,,^,{h\^F^,F,'^ - ^ha^^F.^F^) 

-x'^'ePuefih^ibh^^F^cF/ - \h\bha)pF,,F''') 

-N^t^cuibX\)EdF/y (73) 

Notice that this last expression contains derivatives of the Faraday tensor 
which cannot be replaced by means of the Maxwell equations. These deriva- 
tives enter into the principal part of the evolution equations and destroy the 
hyperbolicity of the evolution equations for the magnetic part of the Weyl 
tensor. In order to deal with this difficulty, in [8] an additional variable, cor- 
responding to the derivative of the Faraday tensor has been introduced. This 
will be discussed in subsection 4.5. 

In order to obtain the contribution of the Klein-Gordon field to the evo- 
lution equation of Bab, we use equation (10) in equation (72). The last two 
terms of equation (72) with Tab replaced by T^^ are a combination of the fields 
(Ta, 0, Aa and the derivative Va^. Accordingly, we introduce the following field 
variable: 

= "^aCk- (74) 

However, the term 

^V^ in^/hP^^^bhafN^) , (75) 

in equation (72), contains first and second derivatives of the Klein-Gordon 
field. In order to recast the term (75) in a more convenient form, we first 
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evaluate ^cT^'^\ Using equation (10) one can write 

(76) 

Important for our purposes is that the last term in the last equation, namely 
4>S/ c4>{t'i^^ + T^^' + r^b*''), contains only (f) and Lpa- Now, consider a generic 
energy momentum tensor 



where C(x) is a constant. A computation then shows that 



VcT^fc' = C(x) ( \{VcXaYb + XaVcYb + V cXbYa + X^V cYa) 

-\gabg'fVcXeYf - ^gahg'fX.VcY^ . (77) 



Applying this last formula to the energy-momentum tensors in (11)-(13) one 
finds that 

'^cT^b^ = C(£r) ((fTaVcO-fe + (TfcVc(T<j) - gabQ^^ <^e^ cC^ f) , (78) 

^cT^b = ((A„V,A + AbVcAa) - ffabs'^^V.A/) , (79) 

^cT^f = C(^) {{(f>a^cA + 0bVc0a) " gabg"^ ^e^ c^^ f) , (80) 

and, moreover, that 

Vc^ir' = Q(^;,VcO-a + AaVcOb + CTa'^cAb + <7bV cAa) 

-^gabg^HAf'^c<Je+<JeVcAf)y (81) 

In order to write this last expression in a form suitable to our purposes we 
introduce the following two auxiliary fields: 

V'ab = VaCTfe, Cab = Va^fe. (82) 



By introducing these new variables we remove all explicit derivatives from the 
term (75). Of course, the price paid by this is that we have to find suitable 
evolution equations for the new auxiliary fields. 
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4.4 Evolution equation for Eab 

The evolution equation for the electric part of the Weyl tensor can be obtained 
by expanding the components P(^a\N\b) ~ ^^abh^^F^Nv of the Friedrich tensor. 
A lengthy computation shows that 

= -^(a|Af|6) ^ 'jl'^abh^^ F'^Nv 

— C-NEab + F>cBd{aeb)''^ " '^O'ce'^'^ (aBb)d ~ '^X{aFb)c ~ '^X'^ (aFb)c 



+habX'''E,, + 2xEab + I F.'^;> - -habh-^FS: ] , (83) 



(a|Af|fc) 2 



where the matter contribution is given by 

{a\N\b) 2 ab"- ^ uNv 

= -tkCn [Tdcha'^'hl 



+KT„,iV''Q/i,,a" - /i(b"a,)^ + \nT^.,N^N-'(^X(ab) - ^X 

+ iAt(^2X(a"^'' - X\aW - Ix^'^habYuv (84) 

In the case of the electromagnetic field the above expression yields explicitly 
that 

77, (cm) it'ab'^ (cm) 

-^(alWIb) 9 ^uNv 



kCn ( h\ah\)FqfFj - ^habFgpF^P ] - ^'^ [ '-^^ ^ ^\aVb) ) E-'F,, 



-nE-F,, {h\aab) - \haba^^ - \nF^dF/ (^2x(a"/ib)" - X^ahf - \x''^Kb 

+ \'^X(ab)Fq,,F'^^ - \nE''x{ab) - IxhabF.pF'^P + \^xhabE\ (85) 

As in the previous subsection we observe the presence of derivatives of the 
Faraday tensor which need to be dealt with by the introduction of a new field 
if one is to preserve the hyperbolicity of the equations. 

Finally, using equation (10) we obtain the expression for the scalar field 
contribution. In this case the first two terms of the general expression (84), 
namely, 

\^Cn {rrhtK) + \ (h^Vb) - h^bV^^ T<->iv", (86) 

clearly contain derivatives of the Klein-Gordon energy-momentum tensor. These 
terms can be rewritten using equation (76) and introducing the auxiliaries 
fields of equations (74) and (82). 
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4.5 Evolution equations for the auxiliary field ipabc 

The analysis of the evolution equations for the eleetric and magnetie parts of 
the Weyl tensor led us to introduce the covariant derivative of the Faraday 
tensor as further field variable. Accordingly, we set 

V'ahc = yaFbc- (87) 

It has been shown in [8] that, applying a covariant derivative to the Maxwell 
equations (61), commuting covariant derivatives and using the definition (87) 
one obtains the following equations for the tensor ipabc 

where 

1 f 

*tpabc = l^ebc ll^aef ■ (89) 

As il^abc = i^a[bc\^ One Can naturally define its electric and magnetic parts 
respect to Na as 

Ead = i^adnN", Ban = N^^a.t- (90) 

Notice that by construction one readily has that £aN°- = BaN'^ = 0. Projecting 
the equations in (88) along the directions longitudinal and transverse to 
one obtains the following set of evolution equations: 

= -2fe"/l/[aV,]7V'' - eabcdhf^'V^B.'^N'') + /l/5ed, (91) 

B.t^f) = -2Se'^/i/[aVfc]A^'' + eabcdhf^'V^Ee^N'') - h/Ved- (92) 
It is noticed that the corresponding constraint equations assume the form 

D'^Bae = -eabcdN'Se'V^N" + N^Ved, (93) 

D^£ea = eabcdN'^Be^'V^N'' ~ N'^S^d, (94) 
where for case of presentation we have set 

^/d = 2F'^''i?[rf|je|6] — V/ J(J, ^fb = -^f'^'^FeuR'^dab- (95) 

4.6 Remarks concerning the unknowns and gauge conditions 

For convenience, we collect the unknowns discussed in the previous sections 
in the following vector variable: 

W ^ {ea^,rf^''^,Eab,Bab.Aa,<P,Ea,Ba,ea,Ba}. (96) 

As already discussed, the matter model under consideration makes further use 
of field Aa and the phase of # = (j)e^^ . Accordingly one has to introduce the 
following auxiliary fields: 

C^a = Va6', ijjab = '^a(^b, Cab = VaA, (97) 
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Recall also, that in equation (74) we have also introduced the 1-form if a = 
Vaft)- Following the discussion in [9] we use the real scalar field (f> to construct a 
frame adapted to the problem. Accordingly, we assume now that (pa is timclike 
and consider a frame satisfying 

(fia = aeo, \a\ = \ip\. (98) 

In fact, in order to construct the adapted frame for this problem, the gauge 
choice has been based on some specific regularity assumptions on the scalar 
field, requiring the timelike vector of the orthonormal frame to follow the 
matter flow. Firstly, it is assumed that (j) G C°°{A4) and, secondly we assume 
'Va4' is and remains timelike. In fact, if the gradient is null or spacelike the 
system evolution breaks down. Hence, one can write Lp"" = aSg. It follows then 
that 

= 0, e'f = 0. (99) 
In terms of components respect to a coordinate basis one finds that 

= ^. (100) 



4.7 Evolution equation for the fields (f> and ipa 

The evolution equation for the real scalar field (f> can be recovered from the 
definition of ipa- Namely, one has that 

dt(b = iff (101) 

An evolution equation for ipa can be found by using equation (16). One has 
then that 

VVa = 0(fTafT"+^'+2AVa). (102) 

In view of the gauge condition in equation (98), equation (102) is, in fact, an 
equation for a and the connection coefficients of the frame. The corresponding 
equations will be discussed in details in Section (4.11). 



4.8 The evolution equation for the vector potential Aa 

The evolution equation for the vector potential Aa can be readily obtained 
from the expression of the Faraday tensor in terms of the curl of the vector 
potential — cf. equation (6). This equation is almost in hyperbolic form. In 
order to close the system one has to introduce a gauge source function — see 
e.g. [5]. This feature is closely related to the gauge freedom built into equation 
(6) and allows to specify freely the divergence of Aa- One has that 

Vo^ = Fob + Vfc^o, (103) 
V°Ao = C{x) - VM„ (104) 

where (^{x) = VqA". The evolution equations for Fab have already been dis- 
cussed in Section 4.2. 
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4.9 The evolution equation for aa 

A key observation for our purposes is that 9, the phase of the complex scalar 
field ^, does not appear explicitly in the evolution equations hitherto discussed, 
appearing only through its derivative. These derivatives are, in turn, replaced 
by (Ta = Va^. The evolution equation for the vector potential aa can be inferred 
from the following integrability condition 

VaCTfc - Vb<Ta = 0, (105) 

as the connection being considered is torsion-free. The structure of this equa- 
tion is similar to that of equation (6). In order to bring it to hyperbolic form 
it has to be complemented with information about the divergence of da pre- 
scribed by a gauge source function. From equation (105) if follows that 

VotTf, = VfeCro, (106) 
VVo = iA(x)-VV„ (107) 

where ip = Vera. The gauge source function ip is not truly independent. In 
fact, as a consequence of the Klein Gordon equation (17) one finds that 

+ C ^ -J {a' + A") eo- (108) 

Thus, one readily sees that ip can be expressed in terms of ( and further field 
variables. 



4.10 The evolution equations for the auxiliary field tpab and (^ab 

In order to motivate the construction of suitable evolution equations for the 
auxiliary fields ipab and (^at, we consider first the commutator of covariant 
derivatives acting on a generic tensor Xa- One has that 

Vc^dXh - VdVcXb = RX,dXa. (109) 
Making use of the definition of Y^c = ^bXc we rewrite equation (109) as 

ycYdb - ^dYcb = R\cdXa, (110) 

for the tensor Ybc- making use of this equation for Xa = <7a, Aa one obtains 

Vclpdb-'^di'cb^ R^cd'^a, (111) 
VcCdb-^dCcb = R%cdAa, (112) 

were the components of the tensor R°'bcd in the right-hand side of equations 

(111) and (112) can be reexpressed, using equation (27), as an algebraic expres- 
sion involving the Weyl tensor and the matter fields — cf. equations (41)-(42). 
Thus, the right-hand side is a function of v only. The antisymmetry with re- 
spect to the covariant derivatives in the left-hand side of equations (111) and 

(112) suggests that they imply symmetric hyperbolic evolution equations if 
information concerning the divergence of fields ipab and Cab is provided. This 
is done in the next sections. 



16 



Danicla Pugliese, Juan A. Valiente Kroon 



4-. 10.1 Evolution equation for the field (^ab 

For convenience, define the tensor ga = V'^Cca- If 5a is explicitly known in 
terms of other field variables, then suitable hyperbolic equations for (^ab are 
given by 

VoCia — ViCoa = R^Oia^-e, (113) 
^\0a+y'Qa=9a{x), (114) 

The required information about can be deduced from the Maxwell equation 

V°Fafc = V^VaAb - VVb^a = ~cj)^CTb + Ab) . (115) 

Commuting covariant derivatives one obtains from the latter that or 

V^Cafc - Vb^Aa ~ R'\bAc - -02(^6 + Ab). (116) 

Finally, noticing that ({x) = V"^a and V^WaAb = V^Cafc = gb one obtains 

9b = VbC + i?"'afc^c - 0^(fT6 + Ab). 

Thus, the diverg cncc V^^ca can be rewritten in terms of the derivatives of the 
gauge source function (which are explicitly known) and an algebraic expression 
of matter fields. 

4- 10.2 Evolution equation for the field tpab 

The analysis of the equation (111) is similar to that of equation (112). As- 
suming that the divergence ha = ^'^tpac is known, one readily obtains the 
symmetric hyperbolic system 

Voi>da -'^di'Oa = R^aOdO-e, (117) 
VoV-Oa + ViV'm = ha. (118) 

Now, exploiting that ipab = V'(a6) > one has that 
Hence, one finds that 

hb = VbV + R''"'ab(Te- 

The derivative Vbip can be computed, in turn, using equation (108). Hence, 
one finds that hb can be expressed in terms of the derivatives of the gauge 
function Q and auxiliary fields. 
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4.11 Evolution equations for frame coefficients and connection coefficients 

In order to construct evolution equations for the tetrad and connection coef- 
ficients we follow closely the procedure indicated in [9]. 

First of all, consider equation (24) with = eo = ai^ and = aeg, 
Lp^ = In other words, we require the timelike vector of the orthonormal 

frame to follow the "matter flow lines Vt/i", — see Section 4.10 and Section 4.6. 
Moreover, we introduce a coordinate system (a;^) by requiring the Lagrange 
condition = c-a ^ — Section 4.1). Taking into account the latter 
and using equation (22) we find that 

a,ef = (V,-r,>^^ + ro°e(;. (ii9) 

The later is the required evolution equation for . 

We turn now the attention to the evolution equation for the connection 
coefficients r^'^f^. The requirement that the frame components arc Fermi 
propagated along eo readily implies that F^^ = — see e.g. [6]. Thus, no 
evolution equations are required for these components. The only components 
requiring evolution equations are ^q^j, ^o'^i ^^'^ j- 

The evolution equation for F^"^ ^ is found using equation (26), where 

dtF^''^ = R^jaq - r^\jF^\ - Ff^\F°j + Fq°^F^\. (120) 

Now, using equation (43) one can express the component of the Riemann 
tensor R''jOq in terms of the matter fields. Using the Einstein equation i?op = 

p(KC 



kTqp, with Top = ri!"' -I- T^^^^ and 



where 

T^^^ = do<j>dp<p ^ 0, F^;^^2aoap, T^^> = 2AoAp, T^^;*) = 4^(0^^(122) 
we have 

i?op = - g'^^FoaFbp + 0^(CTocrp + AfiAp + 2CT(oAp))y (123) 

Accordingly, the component R^p can be expressed purely in terms of undiffer- 
entiated field variables. 

The evolution equations for the coefficients 7})° and 7^^° can be found 
again using equation (26). One finds that 

Now, using the conservation equations and the condition V[aVh]0 = one has 
the symmetry F^ — Fj\. From the Klein Gordon equation for the field (j) we 
can recover the following evolution equation for the variable a: 

d-ax-0(cr^+^^+2crMa) =0. (125) 
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Now, the evolution equation for can be found observing that 

[eo, e,]a = r^\dta - V^a. (126) 

We can use equation (125) for the term dta. From the commutator [e^, eo]0 = 
we obtain the constraint 

V,a - = 0. (127) 

Therefore equation (126) can be rewritten as 

[eo, e,]a = F^^ax - x/(/^o° ") + + + 2gaM,), (128) 

so as to obtain 

- (Vo^o", + v.x) = r,\x - x/r^, + -r,%^ + -v,^, (129) 

■'a acp a acp 

where 

^ = 0(^2+^2^2^'^^,). (130) 

The last term of equation (129) contains Vicr" and ViAa- However, making 
use of the auxiliary variables ipat and (^ab, it can be expressed in terms of 
undifferentiated quantities. Thus, noting that the component i?op of the Ricci 
tensor is 

-^op = ^[ep/^i — GiPpQ + rj^f^lFp^ — i^j p) + -Tj'o-^i p — -^j j-^p'^ol' (131) 
and using equation (131) in equation (123) it follows 

~ \^P^i ~ Si^p + ai-^p^ ~ p) + p ^ j-^p'o\ ~ 

K [-g''''FoaFbp + (/)2 {aoap + AoAp + 2cr(oAp))] . (132) 

Substituting in equation (129), the term (— V^x), where x = —h'^'^T^ch ^^^^ 
noting that —ep^^^Q = Cph'^^ F^^j we have, from equation (132) 

VpX = -e.F;, + r/o(r/, - F^ g + F/,F^ ^ - F, )F^^, + (133) 
K [-g''''FoaFbp + 02 (aoCTp + AoAp + 2ct(o^p))] , (134) 

thus equation (129) reads 

Voro° - epF.^l ^ ^iF,%x - x/ro°. + ^ro",g + ^V,g) (135) 
" {^/o^p'i - ^P^'j^. + ^ [-a^'FoaFM + 02 (aoa, + AoA, + 2a(oA,))] } , 

where we used F^^j = 7^^°. Equation (135) is therefore the evolution equation 
for the connection r'n°-. 
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We concentrate now our analysis on equation (124) for the case of the 
components r^^y Using equation (27), the curvature component i?jOg can be 
spht as 

-RjOg = '^°jOg + 2 {^qj + gjqSj') ■ (136) 

Thus, one finds that 

(Jt^ q j tiql g J — O jOq ^ \^qj + yjq^O ) ~r ^ q ^ i q ^0 q ji 

(137) 

where 

Sab = ^ir' + ^Ir' ' Sab = l^{Tab - lgab9'"'T,d) ■ (138) 

Using equations (4), (9) and (121)-(122), one can readily conclude that that 
the right hand side of equation (137) does not contain derivatives of the field 
unknowns. 

Summarizing, equations (135) and (137) provide the required evolution 
equations for the connection components rjj^ and respectively. 



4.12 Summary of the analysis 

We now summarizes the analysis of the evolution equations carried out in the 
preceding sections. The unknowns of the evolution equations can be conve- 
niently grouped in the vector 

V = (e/, /"o" , rq'-j,rq°j,Eab, Bab, Ea,Ba,£a, Ba, (j), ^a,Cra,a, tpab, Cab)- (139) 

The components of Fa'^c not included in this list are determined by means 
of gauge conditions and symmetries. By construction, the electric and mag- 
netic parts of the Weyl tensor are tracefree. This symmetry is disregarded in 
these considerations and is recovered by imposing it on the initial data. It can 
be shown that if these tensors are initially tracefree, then they will be also 
tracefree for all later times — see e.g. [7]. 

The evolution equations for the independent components of the unknowns 
in (139) and the underlying assumptions in their construction are given as 
follows: 

(i) The evolution equation for the tetrad coefficients, e,;^, is given by equation 
(119) by virtue of the Lagrange condition iV" = = Sq. Equation (119) 
has the same principal part than the corresponding equation in the case 
of an uncharged perfect fluid — see e.g. [4]. It gives rise to a symmetric 
hyperbolic subsystem of equations. 
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(ii) The the evolution equations for the connection coefficients, and /^^ ^ 
are given, respectively, by equations (135) and (137). As a consequence 
of the gauge conditions one has, in addition, that rj^f = 0. Equation (120) 
also takes care of F-'^.. As in the case of the equations frame coefficients, 
equations (135), (137) and (120) have the same principal part as those 
of the analysis given in [4]. Again, one has an hyperbolic subsystem of 
equations. 

(iii) The evolution equations for the electric and magnetic parts of the Weyl 
tensor, Eab and Bab, are given, respectively, by equations (83) and (71). As 
mentioned before, the tracefreeness of these tensors is not used to reduce 
the number of independent components. Thus one has 12 equations for 
equal number of components. Equations with a principal part of the form 
of equations (83) and (71) arc symmetric hyperbolic independently of the 
gauge choice — see e.g. [7]. 

(iv) The evolution equation for the electric and magnetic parts of the Faraday 
tensor, Ea and Bf,, are given, respectively, by equations (67) and (68). As 
in the case of the equations for the electric and magnetic parts of the Weyl 
tensor, the principal part of these equations is known to be hyperbolic 
independently of the gauge — again, see e.g. [7]. 

(v) The evolution equations for the electric and magnetic parts of parts of the 
auxiliary field ipabc (encoding the covariant derivative of the Faraday ten- 
sor), £a and Bb are given, respectively, by equations (91) and (92). These 
equations involve 24 equations for as many unknowns. Their structure is 
analogous to that of equations (67) and (68), except that they contain 
an extra free index. As a consequence, their principal part gives rise to a 
symmetric hyperbolic subsystem. 

(vi) The evolution equation for the vector potential Aa is given by equations 
(103) and (104). The symmetric hyperbolicity of the equations is imple- 
mented by the addition of the (arbitrary) gauge source function (^{x) which 
allows to freely specify the divergence Va^°- 

(vii) The evolution equation for the scalar field (p, is given by equation (101). 

(viii) The evolution equations for the auxiliary field aa = Va^ are given by 
equations (106) and (107) with the understanding that the derivative 
Va(T° — ip{x) arising in the right-hand sides is recovered from the Klein 
Gordon equation in the form given by (108). 

(ix) The equation for the scalar a defined by the conditions a'^ = ip and 
cp'^ = Va(/>V''(/) is given by equation (125). This evolution equation is 
a consequence of the Klein Gordon equation for the field (j) — see equation 
(16). 

(x) The evolution equations for the components of ipab = ^a'fpb are given by 
equations (117) and (118). 

(xi) Finally, the evolution equations for the components of ^ab = Va^ti are 
given by equations (113) and (114). 



Summarizing, equations. (119), (135), (137), (120), (83), (71), (67), (68), 
(91), (92), (103), (104), (106), (107), (125), (117), (118), (113) and (114), 
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provide the desired symmetric hyperbolic evolution system for the Einstein- 
charged scalar field system. This system can be written in a concise form using 
matrix notation. One has that 

A°5ov- A^'9jV = Bv, (140) 

where A'' = A''(a:^,v) are matrix- valued functions of the coordinates and 
the unknowns v, and A" is positive definite at least close to a fiducial initial 
hypersurface. The structure of the system (140) ensures the existence of local 
solutions to the evolution equations. The analysis of whether these solutions 
give rise to a solution of the full Einstein-charged scalar field system requires 
the analysis of the evolution of the constraint equations. This is a computa- 
tionally intensive argument which will be omitted here. A general argument to 
handle the propagation of the constraints without having to resort to lengthy 
computations can be found in [1]. 



5 Conclusion 

In the present article we have revisited the issue of wcll-posedncss initial value 
problem for the evolution equations of the Einstein-Maxwcll-Klcin-Gordon 
system (a self-gravitating charged scalar field). The approach followed makes 
use of the well known 1 -f 3 tetrad formalism by means of which the various 
tensorial quantities and equations arc projected along the direction of the co- 
moving observer and onto the orthogonal subspace. Following [4,7], we require 
the timelikc vector of the orthonormal frame to follow the matter fiow lines 
(Lagrangian gauge) . The adapted frame for this problem has been constructed 
taking cq = 'fia/oi. The gauge condition depends therefore on the scalar field 
evolution, the construction of the appropriate frame for this problem relies 
on suitable regularity conditions for 0. As discussed in Section 4.6 it is as- 
sumed that (f) e C°°{M.) and Va(/> is timelikc. Moreover, we assume the vector 
fields tetrad to be Fermi transported in the direction of U , these conditions 
fix certain components of the connection. 

A key feature of our analysis was the introduction of several auxiliary 
fields, the tensors 'ipabc,Cab,'4'ab corresponding to the covariant derivative of 
the Faraday tensor Fab, the vector potential Aa and the scalar field </>. The 
purpose of introducing these tensors was to ensure the symmetric hyperbolicity 
of the evolution equations for the components of the Weyl tensor. 
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